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Abstract 

Objective: The aim of our study was to investigate the association between the TNF-!x-308G/A polymorphism and 
obstructive sleep apnea syndrome (OSAS). 

Method: The IVledline, Web of Science, ElVlBASE, Chinese National Knowledge Infrastructure (CNKI), and Cochrane Central 
Register of Controlled Trials were searched. Pooled odds ratios (ORs) and 95% confidence Intervals (CIs) were calculated to 
study TNF-a.-308G/A polymorphism and risk of OSAS. 

Result: ^0 case-control studies were Included in our meta-analysis. The results from our study showed that the TNF-a-308G/ 
A polymorphism was significantly associated with risk of OSAS (A vs. G: OR =1.67, 95% CI = 1.43-1.95). In the subgroup 
analysis by ethnicity, the statistical similar results were observed both in European (A vs. G: OR =1.68, 95% CI = 1.35-2.08) 
and Asian population (A vs. G: OR = 2.02, 95% CI = 1.50-2.71). When stratified by age, a significantly Increased risk was 
observed in adult carries A allele compared with G allele (0R= 1.79, 95% Cl = 1.50-2.13), whereas no association was found 
In children (OR =1.09, 95% CI = 0.70-1.69). 

Conclusion: Our study suggested that the TNF-a.- 308C/A polymorphism contributed to the susceptibility to the risk of 
OSAS. Additional well-designed large studies are needed to validate our findings. 
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Introduction 

Obstructive sleep apnea syndrome (OSAS) is a common sleep 
disorder characterized by repetitive partial or complete obstruc- 
tion of the upper respiratory tract during sleep, resulting in apnea 
or hypopnea [1]. Due to obesity and aging population, the OSAS 
has undergone an increasing prevalence all over the world. It was 
reported that more than 5% of the general population has been 
affected [2]. According to the National Sleep Foundation (NSF) 
Sleep in America, there were 1/4 Americans at high risk of 
suffering sleep apnea on the basis of the Berlin Questionnaire [3] . 
OSAS has been reported to be associated with various health 
related consequences, including cardiovascular disease, hyperten- 
sion, stroke, insulin resistance and all-cause mortality [4]. OSAS 
represents a vital public health concern and should be given much 
more attention because of the high prevalence and its enormous 
negative consequences. In consequence, improving our under- 



standing of the pathogenesis of OSAS is essential for the 
development of effective and safe therapies. 

Tumor necrosis factor (TNFl-a, a member of the TNF/TNFR 
cytokine family, is an intercellular communicating molecule 
involved in a wide variety of human diseases. Krueger et al. [5] 
has pointed out that TNF-a is one of the most important 
pleiotropic proinflammatory cytokines involved in sleep regula- 
tion. Raised levels of circulating TNF-ot in patients with OSAS 
have been reported in previous studies [6-8]. The synthesis of 
TNF-a has been suggested to be mostly regulated at the 
transcriptional level [9]. The DNA variations in the promoter 
region of the TNF-a gene may directly influence the transcription 
of the TNF gene. The TNF-a. gene is located within the highly 
polymorphic major histocompatibility complex (MHC) region on 
the short arm of chromosome 6p21.3 [10]. Several polymorphisms 
in the promoter region of the TNF-a gene have been identified. 
Among which, polymorphism at position —308 in the promoter 
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Potoitial relevant studies 
identified by electronic 

databases seaixh update to 
march 31, 2014 (n = 88) 



14 studies for further 
screening 



10 studies finally included 
in the meta-analysis 



74 studies excluded by title or abstract screening: 
Reviews or meta-analysis: n = 4 
Not relevant to TNF-a-308G/A 
polymorphism and OS AS: n = 70 



4 studies excluded by full text assessment: 
Duplicate publications: n = 2 
No interested outcomes: n = 2 



Figure 1. Flow chart of literature search and study selection. 

doi:l 0.1 371/journal.pone.01 061 83.g001 



region, consisting of a guanine (G) to adenine (A) substitution, has 
been reported to be associated with increased production of TNF- 
a levels both in vitro and in vivo [11,12]. There is evidence that A 
allele is associated with increased level of TNF-a in plasma 
compared with G allele [13]. As the TNF-a-308G/A polymor- 
phism is strongly associated with circulating TNF-a concentra- 
tions, it is assumed that it might he closely related to the OSAS 
risk. To investigate a possible association between TNF-a 
polymorphism and risk of OSAS, we conducted a meta-analysis 
from all available relevant studies. 

Materials and Methods 

Literature search strategy and eligibility criteria 

The Medline, Web of Science, EMBASE, Chinese National 
Knowledge Infrastructure (CNKI) and Cochrane Central Register 
of Controlled Trials were searched. A broad search strategy was 
used for optimum sensitivity. Using Medical Subject Headings 
(MeSH) and text words, we adopted the following terms to search 
the databases: ("obstructive sleep apnea-hypopnea syndrome" 
OR "obstructive sleep apnea" OR "sleep apnea" OR "apnea" 
OR "OSAS" OR "OSA") AND ("single nucleotide polymor- 
phisms" OR "SNP" OR "polymorphism" OR "gene variant" 
OR "mutation") AND ("tumor necrosis factor a" OR "tumor 
necrosis factor-a" OR "tumor necrosis factor" OR "TNF-Qf" OR 
"TNF"). The search was restricted to humans. Tides and abstracts 
were screened up to 31 March 2014 were retrieved. Articles were 
screened at the title and abstract phase by two authors (Yanping 
Wu and Chao Cao). 

Inclusion and exclusion criteria 

Inclusion criteria for studies included: (a) evaluation of the 
relationship between TNF-a- 30SG/A polymorphism and OSAS 
susceptibihty; (b) case-control study; (c) validated genotyping 
methods were used; (d) detailed genotype frequencies in cases and 
controls for the calculation. Major reasons for exclusion of studies 
were: (a) no control group, case reports or review articles; (b) 
duplicate data; (c) enrolled patients were using TNF inhibitors. 



Data extraction and quality assessment 

Information was carefully extracted from all eligible studies 
independentiy by two of the authors, according to the inclusion 
criteria listed above. From each study, data on general study 
characteristics was extracted: first author's surname, year of 
publication, ethnic descent of the study jjopulation (European, 
Asian or mixed), genotyping methods, numbers of eligible cases 
and controls, genotype distributions in cases and controls and 
available subgroups. Disagreements were resolved by discussion 
and consensus. 

The Newcasde-Ottawa Scale (NOS), which was developed to 
assess the quality of non-randomised studies, was performed to 
assess the quahty of included studies [14]. The system included 
three broad perspectives: the selection of the study groups; the 
comparability of the groups; and the ascertainment of the 
outcome. And the scores of 0-3, 4-6, and 7-9 indicated low, 
moderate, and high quahty of studies, respectively. 

Statistical analysis 

All analyses were performed on Review Manager version 5.0. 16 
and Stata version 12.0. The Hardy-Weinberg equilibrium (HWE) 
was utilized to compare the observed genotype frequencies with 
expected genotype frequencies in controls. The association 
between risk of OSAS and TNF-a.-308G/A polymorphism was 
estimated for each study using the odds ratio (OR) and 95% 
confidence interval (CI). Between-study heterogeneity was assessed 
with the x^-based Q_ statistical test [15]. We also used the statistic 
of to eflBciendy test for the heterogeneity, with /^<25%, 25- 
75% and >75% to represent low, moderate and high degree of 
inconsistency, respectively [16,17]. When P>0.1 or /^<50%, 
indicating a lack of heterogeneity, and for these analyses, the fixed- 
effect model was used; otherwise, the random-eflfei:t model was 
applied [16-19]. The significance of the combined OR was 
determined by the Z-test, in which _P<0.05 was considered 
significant. Effects were tested in models as follows: allele 
comparison (A vs. G), homozygote comparison (GG vs. AA), 
heterozygote comparison (AG vs. AA), dominant model (AG+GG 
vs. AA), and recessive model (GG vs. AA-I-AG). Stratified analysis 
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OSA 



Control 



Peto Odds Ratio 



Studv or Subarotip 


Events 


Total 


Events 


Total 


Weltiht 


Peto, Fixed, 95% CI 


Almpanidou 2012 


121 


440 


112 


638 


28.4% 


1.80 [1.34,2.41) 


Bhushan 2009 


38 


208 


16 


206 


7.6% 


2.52 [1.42, 4.46) 


Guan 2013 


107 


1062 


20 


324 


13.3% 


1.60 [1.04, 2.46) 


Karkucak 201 2 


18 


138 


9 


84 


3.6% 


1.24 [0.54, 2.84) 


Khalyfa2G11 


47 


276 


48 


302 


12.8% 


1.09 [0.70,1.69) 


Li 2006 


16 


48 


14 


96 


3.4% 


3.09 [1.32, 7.24) 


□ 2013 


18 


310 


12 


200 


4.3% 


0.97 [0.45, 2.05) 


Liu 2006 


39 


152 


8 


84 


5.6% 


2.74 [1.41, 5.33) 


Popko 2008 


29 


204 


18 


154 


6.5% 


1.25 [0.67, 2.32) 


Rilia2005 


58 


206 


68 


380 


14.5% 


1.83 [1.21,2.77) 


Total (95% CI) 




3044 




2468 


100.0% 


1.67 [1.43, 1.95] 


Total events 


491 




325 








Heterogeneity: Chi= = 


13.65, df 


= 9(P = 


0.14); 1= 


= 34% 







Peto Odds Ratio 
Peto, Fixed. 95% CI 



Test for overall effect: Z= 6.40 (P < 0.00001) 



0.01 0.1 1 10 100 
Favours OSAS Favours Control 



Figure 2. Meta-analysis with a fixed model for the ORs of OSAS risk associated with TNF-a-308G/A polymorphism (A vs. G). 

doi:10.1371/journal.pone.0106183.g002 



was performed by ethnicity and age. Sensitivity analysis was 
performed to evaluate the effect of each individual study on the 
pooled ORs. Funnel plots and Egger's linear regression test were 
used to evaluate the potential publication bias [20] . 

Results 

Cliaracteristics of studies 

The search strategy retrieved 88 potentially relevant studies. 
After initial search by the tides and abstracts, 74 studies were 
excluded. Among the excluded studies, 4 were reviews or meta- 
analysis and 70 were not relevant to TNF-a-308G/A polymor- 
phism and OSAS. By further full text screening, 2 studies were 



duplicate publications and 2 had no interested outcomes. Finally, 
10 full-text articles, with a total of 1522 OSAS cases and 1234 
controls, which definitely addressed the association between TNF- 
a.-308G/A polymorphism and OSAS risk were enrolled in our 
study [21-30] (Figure 1). According to the NOS quality assess- 
ment, all included studies achieved a score of moderate to high 
quality (table SI). The genotype distribution in each study was 
extracted and detailed information in one study was referred to 
pubhshed data. Of aU eligible studies for TNF-a-308G/A 
polymorphism, four were performed in European population 
[23,24,27,30], five in Asian descent [21,22,26,28,29] and one in 
mixed populations (European and Asian descents) [25] . In terms of 



OSAS Control Odds Ratio Odds Ratio 

Studv or Subgroup Events Total Events Total Weight Peto. Fixed, 95% CI Peto, Fixed, 95% CI 

4.8.1 European 



Almpanidou 2012 


121 


440 


112 


638 


53.6% 


1.78 [1.33,2.39) 


m 


Karkucak 201 2 


18 


138 


9 


84 


6.3% 


1.25 [0.53, 2.93) 




Popko 2008 


29 


204 


18 


154 


11.6% 


1.25 [0.67, 2.35) 




RIha 2005 


58 


206 


68 


380 


28.5% 


1.80 [1.20,2.69) 




Subtotal (95% CI) 




988 




1256 


100.0% 


1.68 [1.35, 2.08] 


♦ 


Total events 


226 




207 










Heterogeneity: Chl= = 


1.57, df= 


3(P = 


0.67); P = 


0% 









Test for overall effect: Z= 4.73 (P < 0.00001) 



4.8.2 Asian 



Bhushan 2009 


38 208 


16 


206 


22.8% 


2.65 [1.43, 4.93) 




Guan 2013 


107 1062 


20 


324 


35.8% 


1.70 [1.04, 2.79) 




Li 2006 


16 48 


14 


96 


12.8% 


2.93 [1.28,6.69) 




Li 2013 


18 310 


12 


200 


15.4% 


0.97 [0.45, 2.05) —' 




Liu 2006 


39 152 


8 


84 


13.2% 


3.28 [1.45, 7.40) 




Subtotal (95% CI) 


1780 




910 


100.0% 


2.02 [1.50, 2.71] 


♦ 


Total events 


218 


70 










Heterogeneity: Chl= = 


7.03, df= 4 (P = 


0.13); l= = 


43% 








Test for overall effect: Z= 4.65 (P < 0.00001) 











0.01 0.1 1 10 100 

Favours OSAS Favours Control 



Figure 3. Meta-analysis with a fixed model for the ORs of OSAS risk in subgroup analysis by ethnicity (A vs. G). 

dol:1 0.1 371/journal.pone.01 061 83.g003 
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age, nine studies were conducted in adults [21-24,26-30] and one 
in children [25]. Polymerase chain reaction-restriction fragment 
length polymorphism (PCR-RFLP) [22-24,26-28] and TaqMan 
assay [21,25,29,30] were used to analyze the genotyping. The 
genotype distributions of all control samples except one study were 
in HWE [26]. The general characteristics and genotype informa- 
tion of each study were summarized in Table 1 and Table 2. 

Quantitative synthesis 

Our study suggested a significant association of the TNFoc- 
308G/A polymorphism with OSAS risk (Aw?. G: OR= 1.67, 95% 
CI= 1.43-1.95, P<0.00001; homozygote comparison: OR = 2.62, 
95% CI= 1.68-4.09, P<0.0001; dominant model: OR=1.74, 
95% CI= 1.36-2.22, P<0.00001; recessive model: OR = 2.15, 
95% CI= 1.39-3.32, P = 0.0006, Table 3; Figure 2). However, no 
significant association between this polymorphism and OSAS risk 
was observed in heterozygote comparison (OR = 1.38, 95% 
CI = 0.86-2.21, P = 0.18). In the subgroup analysis by ethnicity, 
similar results were observed in European (A vs. G: OR = 1 .68, 
95% CI = 1.35-2.08, P<0.00001) and Asian population (A vs. G: 
OR=2.02, 95% CI= 1.50-2.71, P<0.00001; Figure 3). When 
stratified according to age, significant association between the 
TNFa-308G/A polymorphism and risk of OSAS was observed in 
adults (A vs. G: OR= 1.79, 95% CI = 1.50-2.13, P < 0.00001), 
but not for chUdren (A vs. G: OR= 1.09, 95% CI = 0.70-1.69, 
P=0.71). 

Sensitivity analysis 

Sensitivity analysis was carried out to assess the influence of 
each individual study on the pooled ORs. By removing each study, 
the pooled ORs were not altered significandy, which indicated 
that no individual study significandy affected the pooled results. 
When one HWE-violating study was excluded, the corresponding 
pooled ORs were stiU meaningful, showing the results were robust 
(OR= 1.62, 95% CI= 1.37-1.91, P<0.00001). 

Publication bias 

Begg's test and funnel plot were performed to verify the 
publication bias of the literature. The shapes of the funnel plots did 
not show obvious asymmetry (Figure 4). The asymmetry of the 
funnel plot was further tested by Egger's regression (t = 0.37, 
P= 0.724), which showed no significant publication bias. 

Discussion 

In this meta-analysis, we concluded the existing data on the 
association of TNFa-308G/A polymorphism and OSAS risk. 10 
studies involving 1,522 cases and 1,234 controls met the eligibility 
criteria. Our results demonstrated that the TNF-a-308G/A 
polymorphism was associated with OSAS risk (A vs. G: 
OR= 1.67, 95% CI= 1.43-1.95, P<0.00001). The A aUele, 
compared with the G allele, was more likely to promote the risk 
of OSAS. 

In stratified analyses, we observed that the association between 
TNFa-308G/A polymorphism and OSAS risk remained signifi- 
cant in Europeans (Avs. G: OR= 1.68, 95% CI = 1.35-2.08) and 
Asians (A vs. G: OR =2.02, 95% CI = 1.50-2.71). Genetic 
background might have significant influence on the susceptibility 
of OSAS. However, up to date, no studies were performed in 
Africans, which limited our further analysis in those populations. 
At the same time, it was worth emphasizing that TNFoc-308G/A 
polymorphism was contributed to the increase of OSAS suscep- 
tibility in adults (A vs. G: OR= 1.79, 95% CI = 1.50-2.13), but 
not for children (A tm. G: 0R= 1.09, 95% CI = 0.70-1.69). This 
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Begg'sfunnel plot with pseudo 95% confidence limits 
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Figure 4. Funnel plot for publication bias of the meta-analysis of OSAS risk and TNF-a-308G/Apo\ymor'py\\sm (A vs. G). 

doi:10.1371/journal.pone.0106183.g004 



might reflect the fact that mcreased age modified the associations 
and might helpful to understand the pathogenesis of adults' OSAS. 

The results from our study showed that TNFa.-308G/A 
polymorphism associated with risk of OSAS. Compared with G 
allele, A allele was closely related to increased susceptibility to 
OSAS. The significance of this polymorphism could be explained 
by its possible influence on the expression of TNF-ot protein. It was 
reported that the A allele of TNF-oc gene increased TNF-a 
expression and resulted in higher TNF-a production by means of 
regulating adipose tissue metabolism [31,32]. Moreover, TNF-o; 
played an important role in regialatuig the pro-inflammatory 
cytokines during sleep and closely related to the level of high 
sensitivity C-reactive protein (hsCRP) [33]. Shamsuzzaman et al. 
demonstrated that the plasma C-reactive protein was elevated in 
OSAS patients and the magnitude of CRP elevation was 
remarkably related to the severity of OSAS [34] . 

Compared with the previous meta-analysis performed by Huang 
et al [35], our study had several superiorities. Firsdy, the previous 
study identified less than half the number of studies as compared 
with the present study. More subjects were included in our study to 
provide a more comprehensive result. Secondly, in the previous 
meta-analysis, the authors did not evaluate some clinically 
significant outcomes and clinically relevant subgroups. This is first 
meta-analysis showed TNF-a-308 polymorphism increased suscep- 
tibility to the OSAS in both in European and Asian population. In 
addition, we firstly observed TNF-a-308 polymorphism increased 
the risk of OSAS in Adults, but not for children. Thirdly, no 
significant between-study heterogeneity and publication bias were 
observed in our study. Moreover, we carried out sensitivity analyses 
to evaluate the effect of each individual study on the pooled ORs. 

Some limitations of this meta-analysis should be addressed. 
First, the number of pubhshed studies, especially for African 
population and children, was not large enough. Second, one study 
did not agreed with HWE in our analysis [26]. However, pooled 



ORs were not significant affect when the study was excluded. 
Third, in most studies detail information such as the co- 
morbidities were not available, which limited our further 
performed a risk-adjusted analysis. 

Conclusion 

In summary, this meta-analysis suggested that the A allele of the 
TNF-a-308 polymorphism increased susceptibility to the OSAS. 
These results provided the important role of the TNF-a-308G/A 
polymorphisms in the development of OSAS. Additional well- 
designed large studies were needed to validate our findings. 
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